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Terminal Energy Distribution Of 
Blast Waves from Bursting Spheres 
by 

Andrew A, Adame zyk* 
and 

Roger A, Strehlow 

Aeronautical § Astronautical Engineering Department 
University of Illinois at Urbana-Champaign 


Introduction 

One interesting fundamental property of the blast wave generated by a 
non-ideal source is the fraction of the energy which is initially in the 
source region that ends up in the surroundings after the explosion has oc- 
curred, For ideal sources, i,e,, point source, nuclear or high-explosive 
explosions, virtually all the source energy ends up in the surroundings. 

One can relate this to the fact that these are high-energy density and 
power-density sources. In the case of the ideal burst of a sphere of high- 
pressure, inert gas, this is not true, however, and the energy imparted to 
the surroundings can be considerably less than that initially contained in 
the source region. In this report we will first look at two idealized 
limit cases, and then use actual pressure-volume- time trajectories calcu- 
lated for various bursting spheres to deduce a general relationship for 
the terminal energy distribution of a bursting sphere of ideal gas at high 
pressure. 


* Current address. Ford Motor Co., Research Division, Dearborn, Michigan, 
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Approach : 

For any non-ideal explosion, the energy added at the source, Q, eventually 
partitions itself between the source volume as residual source volume energy. 
Eg, and the surroundings. The quantity of energy added to the surroundings. 

Eg, is determined by the work that the source volume does on the surroundings. 
This work must be evaluated at the contact surface between the source volume 
and the surroundings, and is given by the path integral 



V 

o 


where P is the instantaneous pressure at the contact surface, is the initial 
source volume, and is the final source volume that one finds when the system 
is quiescent, and its pressure has a uniform value , P^, equal to the external 
atmospheric pressure. 

In each case that we consider, we will assume that both the source volume 
and surroundings contain an ideal gas which has a constant heat capacity and 
we will approximate, as appropriate, the actual path integral observed from 
numerical calculation of the blast wave by a suitable simplified path so as to 
be able to analytically or numerically evaluate equation 1. In all cases the 
initial pressure (before energy addition], both in the source volume and in 
the surroundings, will have a constant value of P^. However, density, temper- 
ature, molecular weight and heat capacity of the source gas may in general be 
different than that of the surroundings. With these restrictions we consider 
the following cases: 

I. Constant pressure addition of energy 

II. Constant volume addition of energy with subsequent isentropic expansion 

III. Constant volume addition of energy with the subsequent flow calculated 
using the equations of motion, (i.e., the bursting sphere). 
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Case I - Constant pressure energy addition 

In this case we assume that energy is deposited homogeneously throughout 
the source region, i.e., we assume that the rate of energy addition per unit 
mass within the source volume is the same for all elements of the source vol- 
ume. Additionally, we assume that the rate of energy addition is very slow 
- so slow that the pressure remains constant during the energy addition pro- 
cess. With these assumptions, equation 1 may be integrated to yield the 
equation 


Ec = P. (V, 




( 2 ) 


Furthermore, since the energy is added at constant pressure, the final , 
temperature in the ball is given by the equation 


T 


2 



(3) 


where N is the number of moles of gas in the source volume, C^ is the constant 
pressure heat capacity of the gas in the source volume and the subscripts 1 
and 2 refer to the initial and final gas temperature in the source volume, 
respectively. Furthermore, the energy which remains in the source region 
after Q units of energy have been added is 


E 


B 


NC^ (Tj - Tj) 


C4) 


and the energy in the surroundings is therefore 

Q - Eb . Eg . NR (Tj - Tj) - C5) 


Thus, the fraction of the energy remaining in the source volume after energy 


addition is 
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E 


B 


Q 


^ = L- z 0.714 for Yi - 1.4 


and the fraction transmitted to the surroundings is 


( 6 ) 



0. 286 for Yj = 1*4 


( 7 ] 


Here y ^ is the heat capacity ratio of the gas in the source volume. Notice 
that the ratios given by equations 6 and 7 are independent of the initial 
source size, V^, of the relative energy addition Q/NC^T^ = q and of the 
gamma of the surrounding gas, y^ ^ They are only dependent on the heat 
capacity ratio of the source vulume gas. Also note that a majority of the 
energy remains in the source volume and that no shock wave propagates into 
the surroundings. This is because the assumption that the energy is 
added slowly allows acoustic level signals to adjust the position of the 
surrounding gas as the energy is being added. Furthermore, since the ex- 
ternal work was performed at the lowest pressure available to the system, 
the fraction of energy transferred to the surroundings must be less than that 
for any other case. 


Case II - Constant volume energy addition with subsequent isentropic expansion 

For this case, we assume that energy is added instantaneously and homo- 
geneously throughout the source volume at time t = 0, with the source gas 
expanding isentropically, (i.e., slowly against counter pressure) when t>0. 
This energy can be expressed in terms of the state variables of the gas in the 
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ball using the equation 

(P, - P )V 

Q - NC^ (Tj - T^) . ■ . °i (8, 

The right-hand side of equation (8) is Brode's formula [1] for the energy- 
contained in a bursting sphere. Thus, a sphere bursting at a Pressure 

and initial temperature T 2 is represented exactly by the instantaneous 

addition of the quantity of energy Q to a source volume gas which was 

initially at pressure and temperature T^, Also note that the stored 

energy in a bursting sphere is given by (Pj - 

of the temperature of the gas in the sphere. However, since P^V^ = NRT^ = 
constant for any specific source volume, the mass of gas in the sphere is 
inversely proportional to the initial temperature, T^^. 

We now assume isentropic expansion of the gases in the source volume 
against an equilibrium counter-pressure which just matches the source 
volume pressure at every instant of time. Since the surface pressure is 
always equal to the internal pressure and since the process is reversible, 
we obtain the relationship 



Equation (9) is identical to that recommended by Baker (2) and Brinkley (3) 
for determining the energy of a bursting sphere. We now define a useful 
dimensionless variable 



_ 2 _ 
p V 
0 o 


( 10 ) 


The quantity q represents the amount of energy that is added to the source 
volume relative to that initially present in the source volume. Note that 
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for fixed Q, q can only change if changes, and P> therefore, repre- 
sents a source energy density. Also note that 


Ti Po 


( 11 ) 


and that, as q ^ 0,and as q 0, fixed Q and 

Manipulating equation (9) and applying the definitions which have 
been presented, we obtain 


-i[< 

n 


(1 + q) - (1 + q) 


1/Y,‘ 


1 


( 12 ) 


Furthermore. 
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(13) 


If equations (12) and (13) are added, they show that the total energy, Q, is 
conserved. Equations (12) and (13) also show that, for isentropic expansion, 
the fractional energy released to the surroundings is a function of the orig- 
inal source energy density as well as the value of If we take the limit 

of Eg/Q as q ^ <»^we obtain 





(14) 


which is the correct point source limit. 

However, using L’Hopital’s rule in the limit as q 0,we obtain 


11 . ^ . 11 . 

q-K) ^ q-K) ^ 




Y,-l 

■Yi 


( 15 ) 


which is equal to the constant pressure energy addition limit. Thus, the 
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fractional amount of energy that the source volume can contribute to the 
surroundings varies by about a factor of 4, i.e., from 1 to about .28 for 
a - 1.4 source volume gas. 

Case III - The massless bursting sphere initially at high pressure 

Using Erode’ s formula [1] we find that the energy trapped in a bursting 
sphere is identical to that for instantaneous energy addition as discussed 
in Case II. We also see that the work performed on the surroundings due to 
isentropic expansion, as given by Baker's or Brinkley's formula, is always 
less than the total energy as calculated using Erode 's formula. Furthermore, 
since this idealized expansion is an isentropic reversible process, we would 
expect the work which is performed on the surroundings to be less for a real 
sphere burst. This is because at the instant the sphere releases its contents, 
the contact surface between sphere material and the atmosphere will drop to 
some lower pressure due to shock wave - rarefaction fan flow -matching at the 
contact surface. This occurs because in non steady, inviscid flow the momen- 
tum equation requires that no pressure gradients can exist unless they are 
also accompanied by flow velocity gradients. Thus at the instant of burst of 
the sphere, a rarefaction fan which has both velocity and pressure gradients 
in it propagates from the edge to the center of the sphere. This causes the 
gas In the sphere near the boundary to have a lower pressure than the sphere 
pressure and to simultaneously propagate away from the center. Outside the 
sphere, in the surrounding gas, the pressure is increased by compression and 
a shock wave is formed. The flow of the surrounding gas is also away from 
the center of the sphere. Since the process near the contact surface is 
quasi steady in coordinates that move with the contact surface, the require- 
ment is that pressure and flow velocity be balanced across that surface at 
all times. In the spherical case, the pressure at that surface does drop 
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with time because of flow divergence, but at the very start of the process 
the region of interest can be assumed to be so thin that divergence is not 
important and the usual shock wave - rarefaction fan flow matching require- 
ments can be used to calculate the flow velocity and pressure at the instant 
of burst. Finally, the pressure to be used in equation (1) always remains 
less than the initial kernel pressure, and may even drop below during 
portions of the process. However, eventually it must reach as a final 
state. The exact trajectory can only be calculated using numerical techniques 
at the present time. 

This type of calculation has been performed for a number of cases by 
Strehlow and Ricker^ . The output of these calculations was examined to 
determine the (P,V) behavior with time of the contact surface that separates 
the sphere contents from the surrounding air. Typical (P,V) trajectories are 
shown in figures 1 and 2, As can be seen when comparing the isentropic ex- 
pansion curves to those obtained by integrating the non-steady, gas-dynamic 
conservation equations, as done in Strehlow and Ricker, the slope of 
each curve shows a markedly higher value initially and then gradually approaches 
the isentropic value. This behavior indicates that the effective y of the 
gas--the ratio of specific heats which could be used to fit a polytropic 
pres sure -volume curve to the actual pressure-volume curve- -has a changing 
value during the expansion, starting with one which is radically higher from 
that of the surrounding gas, to one which is approximately equal to that of 
the surrounding gas . Even though the initial bursting conditions change for 
the results presented in figures 1 and 2, the general behavior of the slope 
is consistent. From these results, it can be seen that the expansion process 
cannot be adequately approximated as a single isentrope, as in case II, and 
a more detailed evaluation of the expansion process is necessary. 

To learn more about the (P,V) behavior of the contact surface, these 
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data were replotted using log V/Vo and log P/Po as coordinates. Graphs for 
eight different initial conditions are shown on figures 3 and 4 as solid 
lines. It can be seen that the (P,V) trajectories, even though curved, can 
be quite easily approximated as two straight lines in some cases, and as a 
single straight line in other cases, A further examination of the flow 
associated with the bursting process showed that, for those curves that 
showed two slopes, the break in the slope occurred at approximately the 
time when the head of the strong inward propagating rarefaction wave 
generated at the kernel edge by the burst just reached the contact surface 
after reflection from the kernel center. Those cases for which only one 
slope was observed were low temperature bursts, in which the sonic velocity 
of the kernel gas was very low, and in these cases, the rarefaction fan 
did not reach the contact surface before the pressure at the contact 
surface reached ambient pressure. 

The reason for this (P,V) trajectory with time is not entirely clear 
at the present time. The very low initial slope that was observed is obviously 
due to the fact that the only flow divergence in spherical coordinates is re- 
sponsible for the pressure drop. As will be shown later, the second slope is 
very nearly equal to the isentropic exponent for the gas in the sphere and it 
will be shown later that if the gas inside the contact surface (the sphere gas) 
can be assumed to have a single definable thermodynamic state (which it does 
not) this slope should be the isentropic coefficient, This proved to be very 
useful in the analysis. 

These observations of the numerically-generated pressure-volume curves 
have led to a rather simple technique for determing the fraction of the energy 
transmitted to the surrounding gas, E^/Q, for the case of an actual sphere 
burst. The technique consists of modelling the actual log pressure-log volume 
trace using two straight-line segments when necessary and a single straight- 
line segment for the other cases. 
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The pressure at the kernel edge and the rate of volume change with this 
pressure is what must be modelled to calculate Es/Q, Figure 5 shows that an 
expansion fan propagates to the center of the sphere with time, reflects and 
eventually overwrites the kernel edge. Earlier we stated that for this period 
of time the (P,V) trajectory is such that it has a straight line behavior on 
a log P-log V plot with an effective polytropic exponent of k^ which is con- 
siderably smaller than the ordinary polytropic exponent for the gas in the 
sphere. This is due to the dynamic nature of the expansion process and is 
a path dependent behavior associated with the initial flow following the 
busrt of the sphere. 

When transformed back into the P-V plane, the straight line relation- 


ships (in the log P-log V space) are converted to the expression which is 


indicative of an ideal polytropic gas P 

P 

0 



C16) 


where - k^ is the slope of the curve in log P-log V space. 

Using figure 6 to define the important state points during the expansion 
process, a simple expression for E^/Q is derived. Again, from equation 1, 



(17) 


From figure 6, it can be deduced that the pressure-volume relationship from 
state 5 to state 4 can be expressed as 



Here, is the initial volume of the sphere and P3 is the pressure at the 
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contact surface just after the instant of burst and is calculated from the 
steady state shock-rarefaction fan balance condition. Also, a similar expression 
can be used for the curve--from state 4 to the point in which ambient conditions 
are reached. This expression is written as 
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p 

0 

^ V ^ 

/oj 

^2 
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where 



from equation 


19 


and the fact that 




V 
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t^O ) 


1/^2 



f 

— = 1, we obtain 
o 



(19) 


( 20 ) 


which now relates the volume of kernel when it has expanded to ambient con- 
ditions to the initial contact surface pressure and the kernel volume when 
the contact surface is overridden by the reflected rarefaction wave. 

To determine this intersection point, point 4, the velocity-time re- 
lationship of the rarefaction fan and contact surface must be integrated to 
obtain the position of each of these boundaries as a function of time until 
they intersect. In essence, one can obtain an expression for by 

integrating the wave propagating speed, u a, subject to the condition that 
it propagate from the kernel edge to the origin, reflect and again propagate 
to the expanding kernel edge. 
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( 21 ) 


12 




a . dt + 
4 


(^4 ^ dt 

at 


( 22 ) 


where is the position of the kernel edge; 

is the position of the rarefaction wave; 

R^ is the original position of the kernel edge; 

dRi 

is the velocity of the kernel edge; 

a^ is the average sound speed at the wave during its propagation. 

At the point of intersection, R = R^, and 

e f 


dR, 


’ 

dt 1 

a 4 dt + 

( 34 ) dt = - a^tj + 

4 J 
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a 4 dt 
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h 


If T = t^ - t, where a. • T = 
2 1 4 


a. dt 
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then from equation 23 
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C24) 
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and 



Where N= 1, 2 or 3, depending on whether the flow field is planar, cylindrical 
or spherical. However, since no simple solution exists fcsr this truly non-steady 
flow field, the results of Strehlow and Ricker are used to obtain V^/V^. 

To determine the amount of energy transmitted to the surrounding gas from 
state 3 to state 4, equation 18 is substituted into the energy integral pre- 
sented in equation 17, which yields 


4' o 


ES, = P V 


0 o 


(P3/P0) CV3/V„) 

Cv/v^)’^i 


d (V/V^) 


( 26 ) 


Furthermore, if one assumes a thermodynamically definable state, the energy 
which remains in the source region after Q units of energy have been added is 


Eb 4 = NjCy (T 4 - (27) 

Remembering that the total energy deposited originally in the kernel, Q, 
equals the sum of the energy transmitted to the surroundings plus the energy 
remaining in the kernel, and is expressed as 




(T 2 


Ti) 


( 28 ) 


Using expressions 26, 27, 28 and the ideal gas law, we obtain 
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Integrating equation 26 from state 4 to the final state, and summing to the 
value obtained at state 4, we obtain the energy transmitted to the surrounding 
gas during the entire process. 

(1-k 
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where A = 

Again assuming a thermodynamically definable end state, the energy remain 
in the kernel can be expressed as 




(31) 


and 


N. R 


(32) 


Using the ideal gas law, equation 32 can be expressed as 
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( 33 ) 


and solving for V^/V^, we obtain 
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substitution of equation 30 into equation 34 yields the following expression for 
the final volume. 


V. 


A 

[V 


1 

1 

A ^ A 


t — 1 

.i-S 

1 

CM 

f^fl 

l-k2 ' 


A. 

.''oj 


1-k^- 

• l-k^ 

1-k^ I-IC 2 

V 

oJ 



J 

P 

0 


= (1-v) 

o 


and substituting equation 20 yields 
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using expressions 18 and 29, we obtain 
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This expression, used in conjunction with equations 25, 36 and 38 (as obtained 
from Liepman and Roshko ^^^ ) 
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(38) 


where Y 2 ratio of specific heat of the gas initially in the kernel, 

yields values for k^^, k^, and V^/V^. Then applying equation 30, the energy 
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to the surrounding gas is calculated. 

In these calculations, it can be seen that the u + a field within the 
kernel is of primary importance, since it explicitly determines the fraction 
of total energy available to the shock wave. Furthermore, if the sound speed 
in the kernel is very low, the rarefaction wave does not override the kernel 
edge before the pressure at the edge has been reduced to the pressure of the 
surrounding gas. In this case, there is no abrupt change in slope of the log 
P-log V curve and k^^ » 'k^ for the entire process. For this special case, 
equation 36 reduces to 


A 


1/k, 



(1 - Y) 



(39) 


Using this expression, equation 38 and equation 30 with k^^ = k^, the energy 
transmitted to the surrounding gas is calculated. 

It should be noted that, since state 4 and the final end state were 
assumed to be thermodynamically definable in this model, k 2 is equal to y 
along this segment of the P-V curve. Furthermore, one can observe in 
figure 3 that the slope of this segment of the actual pres sure -volume curves 
does indeed have a value very nearly equal to the y of the surrounding gas. 
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Results and Conclusions 

Using both the single-line and double-line theories to evaluate the 
fraction of energy transmitted to the surrounding gas, values of Eg/Q for 
each of the log P/P^ log curves in figure 3 are plotted in figure 7 

versus the initial pressure ratio in the sphere. The smooth curves were 
calculated using single line theory for " O.S, 50, and 5000 and using 

double line theory for T^/T^ = 50. These curves are compared to the 
theoretical values for case 1 (constant pressure addition of energy, 

E /Q = 0.286) and for case 2 (isentropic expansion against counter pressure). 

(4) 

They are also compared to the results of Strehlow and Ricker, It should 
be noted here that P 2 /P^-l = q> the energy density of the sphere gas. 

Notice the large scatter in the results of the Strehlow and Ricker 
calculations. Some of this scatter, particularly the lower values of E^/Q 
obtained from some runs is due to the fact that for these cases, the cal- 
culation was not carried out long enough to obtain terminal values. It 
can be noted on figures 3 and 4 that some calculations were terminated when 
the interface pressure was well above atmospheric. The remaing scatter is 
probably due to the late time accumulation of numerical errors. 

The =0.5 curve on figure 7 should represent the actual terminal 

energy distribution, because for these initial conditions the behavior is 
well approximated by a single line curve. The cases on figure 

7 should not agree with the single line calculation because the reflected 
rarefaction wave will reach the contact surface well before the contact 
surface pressure drops to P^. Thus it is not surprising that the Strehlow 
and Ricker values show better agreement to the double-line theory than to 
the single-line theory. 

Finally, it should be noted from figure 7 that the values of E^/Q for 
a bursting sphere, as predicted by this relatively simple theoretical approach. 
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all lie between the Case I and Case II limit values. Furthermore, at fixed 
sphere energy density, q = infinity 

the E^/Q value ranges from the constant pressure limit to the isentropic 
limit. This is as it should be because if the temperature in the sphere 
is very low, the velocity of sound approaches zero and the contact pressure 
approaches P^, On the other hand, if the temperature in the sphere approaches 
infinity, the velocity of sound becomes very large and the expansion is as 
an isentropic expansion at the contact surface because the sphere pressure 
is spatially uniform at all times during burst. Thus we have shown that 
the E^/Q value for all bursting sphere cases lie between the constant pres- 
sure addition case, in which no shock is transmitted to the surroundings 
and the isentropic case in which the maximum work is being performed on 
the surroundings. 
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Figure !• Pressure-volume behavior at contact surface for a typical bursting sphere 
(Hot sphere; P^/P^ = 100, T^/T^ = 50). 
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Figure 2. Pressure volume plot at the contact surface for bursting sphere, 
(cold sphere; P^/P^ = 100, T^/T^ = 0.5). 
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Figure 3. Plot of pres sure -volume at the edge of an expanding explosive 
kernel (Hot spheres; = 50), 





Figure 5, Schematic representation of the path the expansion wave 
travels before encountering the kernel edge. 
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Figure 7. Summary of results. Terminal values of E /Q versus P^/P for different bursting spheres. 
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